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ABSTRACT 


Several  synoptic  problems,  previously  solved  by  graphical  overlays,  can 
be  solved  with  a  programmable  calculator.  The  HP-41CV  hand-held  calculator  is 
used  to  determine  surface  winds,  the  probability  of  average  rainfall  for  a 
regional  area  (Monterey,  California),  contrail  formation,  and  the  lifting 
condensation  level  from  surface  data. 


I. 


Introduction 


This  report  describes  four  programs  that  use  the  HP-41CV  calculator’s 
programming  capability  in  dealing  with  interesting  meteorological  problems, 
including  determining  surface  winds,  the  probability  of  average  rainfall  for  a 
regional  area  (Monterey,  California),  contrail  formation  and  the  lifting 
condensation  level  from  surface  data. 

The  programs  were  prepared  to  use  the  hand-held  calculator  in  solving 
synoptic  problems.  The  first  problem  is  the  estimation  of  surface  wind  at  a 
point  from  a  regional  sea-level  pressure  analysis  or  forecast.  This  problem 
arose  in  the  synoptic  laboratory  where  wind  forecasts  were  prepared  with 
numerical  weather  prediction  guidance  (facsimile  sea-level  pressure  maps). 
The  immediate  solution  is  a  geostrophic  wind  scale  applied  to  the  facsimile 
chart.  The  HP  program  provides  an  alternative  solution.  The  program 
calculates  the  geostrophic  and  gradient  wind,  as  well  as  the  estimated  surface 
wind,  using  Ekman  boundary  layer  assumptions,  by  simply  entering  latitudes, 
longitudes  and  pressure  readings  from  a  sea-level  pressure  chart. 

The  second  problem  involves  preparing  a  statistical  precipitation 
forecast  for  a  set  of  graphical  tables.  Using  a  tutorial  program  in  a  hand¬ 
held  calculator,  the  need  for  the  graphical  tables  is  eliminated.  The  program 
determines  the  probability  of  average  rainfall  over  the  Monterey,  California 
area  for  the  period  0800  PST  forecast  day  to  0800  PST  of  the  following  day, 
for  the  wet-season  months,  November  through  April.  This  procedure  is  based  on 
the  Renard  (1972)  graphical  objective  technique  for  forecasting  24-h  rainfall 
at  Monterey,  California.  From  the  data  entries  of  500  mb  heights  used  to 
calculate  geostrophic  relative  vorticity  at  Monterey  and  a  point  8°  latitude 
upstream  from  Monterey,  as  well  as  sea-level  pressure  data,  the  program  uses 
a  series  of  pre-established  data  tables  to  arrive  at  a  probability  of  24-h 
rainfall  for  Monterey,  California. 
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The  third  problem  involves  forecasting  contrails  from  upper-level 
sounding  data.  It  was  desired  to  have  the  HP  make  the  decision  from  input 
sounding  data  as  an  alternative  to  plotting  the  entire  sounding  and  using  a 
graphical  overlay.  This  program  determines  contrail  formation  at  a  given 
pressure  level  based  on  the  curves  found  in  the  Contrails  Forecasting  Manual 
of  Chief  of  Naval  Operations  (1964).  A  least  squares  approximation  was  made  to 
the  curves  involving  an  expression  for  critical  temperature  which  is  a 
function  of  temperature,  pressure  and  relative  humidity.  If  the  relative 
humidity  is  estimated  (dew-point  temperature  is  unknown),  there  is  a  series  of 
criteria  to  determine  whether  or  not  contrails  will  occur. 

The  final  problem  is  to  determine  the  lifting  condensation  level  (LCL) 
without  using  a  thermodynamic  diagram.  Inputing  values  of  surface 
temperature,  surface  pressure  and  specific  humidity  to  the  HP  results  in  a 
lifting  condensation  level.  The  LCL  is  the  height  at  which  a  parcel,  when 
lifted  dry-adiabatically,  becomes  saturated.  The  program  uses  the  Clausius- 
Clapeyron  equation  to  arrive  at  a  final  expression  for  the  LCL  that  is 
determined  by  values  of  surface  temperature,  surface  pressure  and  specific 
humidity. 


II.  Computation  of  Surface  Winds 
A.  Objectives 


This  program  will  give  the  forecaster,  either  ship-board  or  at  a  regional 
center,  an  easy  and  accurate  method  of  computing  the  geostrophic,  gradient 
and  estimated  surface  wind  from  a  MSL  pressure  analysis  or  model  .forecast. 
This  method  is  more  complete  than  a  graphical  geostrophic  wind  scale  overlay. 

B.  Principles 

The  geostrophic  wind  speed  and  direction  are  obtained  from  the  east-west 
(u«)  and  north-south  (vg)  components. 


ug 


_1  3P 
fp  3y 


|vg|  =  Ju2g 


V^g 


Vg 


1_  3p 
fp  3.x 


tan  0  g  =  -- 

Ug 


where  ug  and  vg  are  geostrophic  wind  components, 

f  =  2  SI  sin  4>  ;  is  latitude  of  forecast  point, 
p  =  1.3  kg/m3;  typical  sea  level  density  value, 

3x, 3y  =  1°  lat  =  111  km 


Figure  1 

Geostrophic  wind  dTg )  computed  by 
the  finite  differencing  scheme 
with  1°  latitude  grid  spacing. 


The  user  enters  the  latitude  and  longitude  of  the  point  (Northern 
Hemisphere)  where  a  wind  is  desired.  Next,  the  user  enters  four  pressures 
that  are  1°  latitude  north,  south,  east  and  west  of  the  center  point  used  in 
the  finite  differencing  technique: 


3p  _  _p(2i_p(D  .  ap  _  _p2):_p© 

3y  ~  ~222~km  ’  3x  '  222  km 


(See  Figure  1 .  ) 


The  gradient  wind  speed  is  obtained  from  the  geostrophic  wind  speed  and 
the  radius  of  curvature  (R) . 

Iv  i  _  _ [V*] _  where  R  >  0  cyclonic 

*r  .5  +  y/25  +  Vg/fR  R  <  0  anticyclonic 

The  user  enters  a  radius  of  curvature  (in  kilometers)  which  is  representative 

of  the  isobaric  curvature  in  the  forecast  area.  (See  Figure  1.) 

In  a  cyclonic  situation  (R  >  0),  the  magnitude  of  the  gradient  wind  is 

less  than  the  magnitude  of  the  geostrophic  wind  (subgeostrophic) ;  if  R  <  0 
(anticyclonic),  the  gradient  wind  speed  is  greater  than  the  geostrophic 
(supergeostrophic) .  If  the  curvature  is  zero  (R  =  » )  the  magnitude  of  the 
gradient  and  geostrophic  winds  is  the  same.  Regardless  of  the  curvature,  the 
gradient  wind  is  from  the  same  direction  as  the  geostrophic  wind. 

The  estimated  surface  wind  speed  and  direction  are  obtained  from  the 
gradient  speed  and  direction  using  Ekman  layer  assumptions.  Over  the  land, 
the  gradient  wind  speed  is  multiplied  by  a  factor  of  0.6  which  is  an 
approximate  amount  that  the  speed  is  reduced  by  friction.  The  resulting  speed 
is  an  estimated  surface  wind  speed.  The  actual  direction  is  rotated  counter¬ 
clockwise  20°  from  the  original  direction.  Over  the  sea,  the  gradient  wind 
speed  is  multiplied  by  .81  and  direction  is  rotated  by  10°.  The  oceanic  and 
land  frictional  effects  follow  the  observational  summary  of  Sheppard  (1969) 
and  the  discussion  of  Petterssen  (1956)  and  Haltiner  and  Martin  <1957). 
Sheppard  (1969)  and  Mendenhall  (1967)  indicate  10°  is  a  reasonable  estimate  of 
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frictional  veering  over  the  ocean  from  weather  ship  data.  Over  land,  the 
frictional  effects  are  naturally  heavily  dependent  upon  surface  roughness  so 
the  20°  veering  is  only  a  general  estimate. 

C.  User’s  Guide 

1)  Make  sure  the  calculator  and  printer  (NOR  mode)  are  on. 

2)  |  XEQ  alpha  SIZE  alpha  I 

will  prompt  SIZE  _ 

3)  Type  100. 

4)  Enter  magnetic  cards  as  labelled  1-9  - 

5)  After  the  program  (WIND)  is  "WORKING"  and  entered,  it  is  ready  to  be 
executed. 

6)  |  XEQ  |  (alpha  |  WIND  (alpha 

7)  The  prompt  CTR  LAT?  will  appear. 

8)  Press  the  |  R/S  |  button  to  continue. 

(Do  this  each  time  you  answer  a  prompt.) 

9)  The  prompt  CTR  LON?  will  appear. 

Enter  the  longitude  of  the  center  forecast  point. 

10)  Press  |  R/S  | 

11)  Answer  the  prompts  for  the  four  pressures  (mb). 

12)  Press  R/S  .  (Remember  to  do  this  each  time.) 

13)  The  geostrophic  wind  speed  (m/s)  and  direction  will  be  displayed. 

14)  Enter  a  radius  of  curvature  (m)  after  prompt  RADIUS? 

15)  The  gradient  wind  speed  will  be  displayed. 

16)  The  prompt  LAND  or  SEA?  will  appear. 

If  the  forecast  region  is  over  land,  enter  "LAND";  otherwise,  enter 


D.  Examples 

To  illustrate  the  use  of  the  program  the  surface  wind  will  be  estimated 
for  points  A  and  B  in  Figure  2.  The  observed  surface  wind  speed  and  direction 
for  forecast  point  A  at  0000Z  on  26  November  1981  is  10  m/s  from  315° .  For 
point  B,  the  surface  wind  speed  is  20  m/s  from  295°. 

Execution  of  the  "WIND”  program  on  the  HP-41CV  calculator  results  in  a 
surface  wind  speed  of  14  m/s  from  a  northwest  direction  of  319°  for  forecast 
point  A.  Entering  the  latitude/longitude  position  (43N  130  W)  along  with  four 
pressures  1°  latitude  north,  south,  east  and  west  (1007  mb,  1012  mb,  1009  mb, 
1012mb  respectively)  results  in  a  geostrophic  wind  of  20  m/s  from  329° .  A 
radius  of  curvature  of  1,100  km  allows  a  gradient  wind  to  be  computed  (18 
m/s).  Finally,  since  forecast  point  A  is  an  ocean  location,  the  oceanic 
Ekman  layer  boundary  assumptions  are  applied.  The  resulting  estimated  surface 
wind  is  14  m/s  from  319° .  This  estimated  wind  is  in  agreement  with  the 
observed  real  wind  speed  (10  m/s)  and  direction  (315°). 

For  forecast  point  B  (46N  135W)  the  four  pressures  (1008  mb,  1011  mb, 
1006.5  mb  and  1013  mb)  result  in  a  geostrophic  wind  of  24  m/s  from  295°. 
From  the  radius  of  curvature  (1,300  km),  a  gradient  wind  of  21  m./s  is 
calculated.  The  final  estimated  wind  speed  for  this  ocean  location  is  17  m/s 
from  285°  which  compares  with  the  observed  real  wind  speed  (20  m/s)  and 
direction  (295°). 


Figure  2 


Surface  conditions  observed  at  0000Z  26  November,  1981  for 

illustration  of  surface  wind  program. 
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Surface  Winds  Program. 
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37  ENTER! 

137  180.0 

137  GTO  -ACT- 

38  STO  88 

88  RCL  18 

138  STO  18 

183*L8L  ‘ANTI 

39  RCL  87 

39  X>Y’ 

139  GTO  -GRAD* 

139  RCL  11 

40  ENTER! 

98  GTO  -TWO- 

I40*LBL  -NOR- 

198  RCL  25 

41  RCL  08 

91  1/X 

141  360.8 

191  / 

42  / 

92  RCL  89 

H2  STO  13 

192  RCL  03 

43  STO  89 

93  « 

143  GTO  -GRAD- 

193  / 

44  RCL  85 

94  ATAN 

144*LBL  -chks- 

194  STO  29 

45  ENTER? 

95  STO  18 

145  0 

195  RCL  20 

46  RCL  06 

96  GTO  ’GRAD* 

146  ENTER! 

196  ENTEP! 

47  - 

97*L8L  -TWO- 

147  RCL  09 

197  .25 

48  108 

98  1/X 

148  X<=Y’ 

198  + 

49  * 

99  RCL  09 

149  GTO  -EAST- 

199  STO  28 

50  2.22  E05  1 

100  * 

150  270.9 

298  9 
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201  ENTEPt 
292  ►'CL  28 
282  J>V'> 

204  GTO  ‘OK- 

205  -0AB  FAB’ 

290  AVI  EH 
207  GTO  “END” 
20S*Lfl  ■?<" 

20?  SORT 

210  .5 

211  + 

212  STO  29 
212  PCL  it 

214  RCl  2? 

215  / 

21c  S™  p 
217*LEl  -QCT- 
21?  “GRAB.  SPEED ; ’ 
21°  PVIEU 
220  PEE 
*221  PCL  27 

222  VIEW  £7 

223  PSE 

224  -LhHB  OR  SEP’- 

225  PON 

226  PROMPT 

227  POFF 

228  PSTO  V 
22?  -LPNB- 

230  PSTO  INK  X 

231  X*Y? 

232  GTO  -LAMB*  % , 

233  PCL  27 

234  EHTERt 

235  .81 
230  * 

237  STO  32 

238  -SFC  SPEED-" 


.  23?  PVIEW 

240  PSE 

241  PCL  32 

242  UIFH  32 
24?  PCL  18 

244  EHTERt 

245  10 
240  - 

247  STO  33 

248  *SFC  BIR  FROM 

249  AVI  EH 

250  CSE 
25S  PCL  33 

252  VIEW  3? 

253  GTO  -EHP- 
2544LBL  “LPMB- 

255  PCL  27 

256  EHTERt 

257  .6 

258  * 

259  ?T0  30 

260  ’SFC  SPEEIi  * 

261  AVI EH 

262  PSE 
26?  PCL  30 

264  VIEW  30 

265  F:CL*18 

266  EHTERt 

267  ?6 

268  - 

26?  STO  31 
276  -SFC  DIR  FROM' 

271  AVI EH 

272  PSE 
273-  PCL  31 
274  yIEH  3! 

2754LBL  -END- 

276  -END- 

277  .END. 
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III.  Determining  the  Probability  of  Precipitation  for  Monterey,  California 

A.  Objective 

To  determine  the  probability  of  24-h  rainfall  for  a  regional  area 
(Monterey,  California)  for  the  period  0800  PST  forecast  day  to  0800  PST  of 
the  following  day,  for  the  wet-season  months,  November  through  April.  This 
procedure  is  based  on  Renard's  (1972)  honogram  method  and  statistics. 

B.  Principles 

General  Forecast  Procedure.  The  probability  of  24-h  rainfall  is 
determined  by  values  of  geostrophic  relative  vorticity  and  sea-level  pressures 
at  different  locations. 

A  first  forecast  parameter  is  obtained  by  taking  the  sea-level  pressure 
at  Eureka,  CA  (Peureita)  and  the  difference  between  the  500  mb 
geostrophic  relative  vorticity  at  Monterey  (  Cm)  and  a  point  8°  latitude 
upstream  from  Monterey  (  Ca).  Enter  these  values  into  a  pre-established 
statistical  table,  which  is  stored  on  the  calculator,  to  arrive  at  a  forecast 
variable,  Yi .  (See  Figure  4.) 

A  second  forecast  parameter,  Ya ,  is  the  value  that  results  from  taking 
the  sea-level  pressure  difference  between  Monterey  and  Eureka  (APm- Eureka) 
and  Monterey  and  Las  Vegas  (APm-lv).  (See  Figure  5).  These  two  parameters, 
Yi  and  Yz ,  are  entered  in  a  third  table  (Figure  6)  which  results  in  a  final 
value  that  corresponds  to  a  probability  of  average  rainfall  for  the  Monterey 
Peninsula. 

500  mb  Geostrophic  Relative  Vorticity.  To  calculate  the  500  mb 
geostrophic  relative  vorticity,  five  500  mb  heights  are  needed  to  use  finite 
differencing  technique  to  solve  (?m  =  g/f  V2  Z)  where  m  is  the  geostrophic 
relative  vorticity  at  Monterey,  V2Z  is  the  Laplacian  of ‘the  500  mb  height. 
For  example,  in  Figure  3  at  Monterey  (Zm  =  558  dm)  the  four  500  mb  heights 
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are  552  dm  (Zni  500  ob  height  to  the  north),  566  dm  i  Z>  > ,  560  dm  i Z*  ; ,  and 
563  dm  (Zw)  with  a  grid  spacing  of  6°  latitude.  This  results  in  a  geostrophic 
relative  vorticity  of  1.8  x  10-5  sec-1. 


Figure  3 

Computation  of  geostrophic  relative 
vorticity  for  Monterey  '5  m)  and  a 
point  8°  latitude  upstream  (Ca). 


The  500  mb  geostrophic  relative  vorticity  is  computed  for  a  point  8° 
latitude  upstream  (8°  corresponds  to  a  typical  daily  progression  of  short- 
waves  in  the  westerlies  of  22  kt  during  the  winter  season)  from  Monterey  in 
the  same  way.  The  difference  between  the  two  vorticities  (Ca®  -Cm)  is  used 
as  a  forecast  factor. 

Sea- Level  Pressure.  Yi ,  a  first  forecast  parameter,  is  obtained  from 
the  computed  vorticity  difference  and  a  sea-level  pressure  at  Eureka, 
California  (Figure  4).  Y2 ,  a  second  forecast  parameter,  is  obtained  from 
the  difference  between  the  sea-level  pressure  at  Monterey  minus  the  sea— level 
pressure  at  Eureka  (APM-surena)  and  the  sea-level  pressure  at  Monterey  minus 
the  sea-level  pressure  at  Las  Vegas  (APm-lv)  (Figure  5). 

Probability  of  Rainfall.  From  the  two  forecast  parameters,  Yi  and  Hz  ,  a 
value  Y3  is  obtained  which  corresponds  to  a  probability  of  rainfall  for  the 
Monterey,  CA  area  for  the  period  0800  PST  forecast  day  to  0800  PST  of  the 
following  day,  for  the  wet-season  months,  November  through  April.  (Figure  6.) 
The  corresponding  values  between  Y3  and  a  probability  of  average  rainfall  are 
on  data  cards  that  are  entered  during  the  program  execution. 
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Eureka 


Figure  4 


f(Peur*ka,  Cs  -Cm):  first  forecast 
factor  of  the  raincaster. 
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C.  User’s  Guide 

1)  Make  sure  the  printer  is  attached  and  in  the  NOR  mode.  Both  printer  and 
calculator  should  be  on  and  the  calculator’s  memory  cleared. 


will  prompt  SIZE _ _ 

3)  Type  035. 

4)  Load  the  program  cards  (1-18).  Continue  after  the  prompt  "WORKING." 

5)  1  XEQ  1  | alpha  j  RAIN  |  alpha 

6)  Enter  500  mb  heights  (dm)  for  each  prompt.  (For  example,  ZN:  8°  means 
height  at  north  position  at  a  point  8°  latitude  upstream  from  Monterey. ) 

7)  Hit  R/'S  after  each  entry. 

8)  Enter  pressure  values  for  Monterey,  Las  Vegas,  and  Eureka. 

9)  Hit  R/S  after  each  entry. 

10)  A  data  card  prompt  will  occur  after  all  the  information  has  been  entered. 
For  the  first  time,  use  data  set  #1;  the  second  time,  use  data  set  #2, 
etc. 

11)  The  final  number  that  appears  after  all  data  card  insertion  instructions 
is  the  probability  of  average  rainfall  for  the  Monterey  Peninsula. 

NOTE:  If  the  program  ends  with  "OUT  OF  RANGE",  then  the  values  used  are 
either  too  small  or  too  large.  Try  another  set  of  values  that 
will  fit  in  the  stored  data  sets.  ' 


0.  Examples 

Example  A  illustrates  the  use  of  the  program  "RAIN"  in  determining  the 
probability  of  precipitation  at  Monterey,  CA.  At  1200Z  25  November  1981  for  a 
point  8°  latitude  upstream  from  Monterey,  the  500  mb  heights  (in  decameters) 
6°  latitude  to  the  north,  east,  south,  west  and  at  the  center  point  are  541.0, 
543.0,  562.5,  563.0  and  555.0  respectively.  (See  Figure  7.)  Around  Monterey, 
the  500  mb  heights  are  550.0  dm  6°  latitude  to  the  north,  555.0  dm  at  the 
eastern  location,  573.0  dm  to  the  south,  576.0  dm  to  the  west,  and  567.0  dm  at 
Monterey.  These  heights  are  used  to  obtain  the  500  mb  relative  geostrophic 
vorticity  values  used  in  generating  forecast  variables.  The  1200Z  sea-level 
pressures  at  Monterey,  Las  Vegas  and  Eureka  (1018.8  mb,  1006.2  mb,  1017.3  mb) 
are  also  required  to  compute  forecast  variables.  (See  Figure  8a.)  These 
forecast  variables  result  in  a  probability  of  precipitation  of  90%.  Rain  did 
verify  during  the  24-h  period,  as  noted  from  the  1200Z  26  November  1981 
observation  (Figure  8b.) 


Figure  8a 

Surface  observations 
for  1200Z  25  November  1981 


Figure  8b 

Surface  observations 
for  1200Z  26  November  1981 


Table  II 

Example  of  the  Probability  of  Precipitation  for  Monterey,  California 
Program. 


ENTER 


ENTER  ZE  S 


fSTCp  jc  0 


EMTE5  :w  ? 


ENTER  z*  - 


ENTER  >  « 


ENTER  ZE  * 


ENTER  Z3  n 


ENTER  ZN  « 


ENTER  N 


ENTER  PRES  NONT 
1.018.8 

Enter  PRES  lvgs 
1 . 808 . 2000 
ENTEP  PPES  ER^Q 
1  0I\ 3006 
vj=  |.at.M000 
X2=-3?.000P 
X3=  Z.0000 

X4=  14.  (WPP 

INSERT  CP  I'  13 

Yl =-20. PPPP 

INEER*  CPP  ? 

Y2=  2E.0000 
00.0000 

Probability  of 
Precipitation  is  90%. 
Data  verified  from 
observed  case  at 
1200Z  26  November  1981. 


Program 


Table  III 

Probability  of  24-h  Precipitation  for  Monterey,  California  Program. 
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240  GTO  24 
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42  PROMPT 
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242  28.035 
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143  STO  83 
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144  RCL  01 
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45  STO  01 
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342  RCL  12 

392  X=Y? 

442  RCL  06 

492  GTO  99 

293  20.035 

343  GTO  24 

393  GTO  65 

443  1 

493  SF  02 

294  -INSERT  CRD  10- 

344*LBL  24 

394  GTO  26 

444  - 

494  RCL  06 

295  AVIEH 

345  RCL  82 

395*LBl  25 

445  STO  06 

495  1 

296  RDTAX 

346  RCL  08 

396  3 

446  GTO  31 

496  ST0  0° 

297  XROH  28.33 

347  X=V? 

39?  - 

447*LBL  31 

497  XEQ  00 

298*L8L  06 

348  GTO  50 

398  PTN 

448  FS’  08 

49g  X’Y’ 

299  GTO  44 

349  XEQ  25 

399*LBL  26 

449  GTO  32 

499  GTO  66 

308  GTO  24 

350  X=Y? 

408  XEQ  -ROUND- 

458  FS’  01 

590  XEQ  00 

aw 


mi 


'  V  -.* 

■  *  '\<o 


531  X>Y? 

551  48 

681  X=V’ 

502  GTO  67 

552  X>YJ? 

682  GTO  61 

503  XEQ  00 

553  GTO  83 

683  XEQ  00 

504  X>Y’ 

554  GTO  84 

684  X=Y? 

505  GTO  G? 

555  RTN 

685  GTO  62 

586  XEQ  80 

556*LBL  99 

606  GTO  63 

587  X>Y? 

557  PCX 

687  GTO  26 

508  GTO  09 

558  Q8V 

688*L8L  35 

509  XEQ  00 

559  CF  01 

689  RCL  86 

510  X>Y'? 

568  CF  82 

618  *v3=‘ 

511  GTO  78 

56!  RTN 

611  QCQ 

512  XEQ  80 

5624181  33 

612  PCX 

513  X>Y? 

563  -2 

613  PSV 

514  GTO  71 

564  STO  09 

614  GTO  23 

515  XEQ  80 

565  RCL  04 

615*LBL  44 

510  X>Y’ 

566  RCL  88 

616  -5 

517  GTO  72 

56?  X=Y'> 

617  STO  09 

518  XEQ  00 

563  GTO  58 

618  RCL  12 

519  JOY? 

569  XEQ  00 

619  RCL  08 

528  GTO  73 

578  X=Y“> 

620  X=Y? 

52!  XEQ  00 

571  GTO  51 

621  GTO  58 

522  \or> 

572  XEQ  00 

622  XEQ  00 

523  GTO  74 

573  X=Y? 

623  X=Y? 

524  XEQ  80 

574  GTO  52 

624  GTO  51 

525  X>Y? 

575  XEQ  80 

625  XEQ  00 

526  GTO  75 

576  XsY? 

626  X=Y? 

527  XEQ  00 

577  GTO  53 

627  GTO  52 

523  x>r? 

573  XEQ  08 

628  m  80 

529  GTO  76 

579  X*Y? 

629  X*Y? 

538  XEQ  00 

580  GTO  54 

639  GTO  53 

53!  XEQ  08 

581  XEQ  08 

631  XEQ  00 

532  t)T' 

532  X=Y? 

632  X=Y? 

533  GTO  7? 

533  GTO  55 

633  GTO  54 

534  XEQ  00 

584  XEQ  00 

634  XEQ  00 

535  XEQ  00 

585  X=Y? 

635  X=Y? 

536  X>f> 

536  RCL  26 

636  GTO  55 

537  GTO  74 

537  GTO  56 

637  XEQ  00 

538  XEQ  00 

533  XEQ  00 

638  X=Y‘J 

539  XEQ  00 

539  X=Y? 

639  GTO  56 

548  X>Y/> 

598  GTO  57 

648  XEQ  80 

541  GTO  38 

591  XEQ  00 

641  X=Y? 

542  XEQ  00 

592  X=Y? 

642  GTO  57 

543  XEQ  80 

593  GTO  58 

643  XEQ  00 

544  X>Y? 

594  XEQ  80 

644  X=Y? 

545  GTO  8! 

595  X=Y^ 

645  GTO  58 

546  CLX 

596  GTO  59 

646  XEQ  00 

547  29 

597  XEQ  80 

647  X=Y’ 

548  JOY? 

598  X=Y’ 

648  GTO  59 

549  GTO  82 

599  GTO  68 

649  XEQ  80 

558  CLX 

688  XEQ  88 

658  X=Y’ 

651  GTO  60 

78!  GTO  26 

?51  75 

652  XEQ  00 

782*l8L  62 

752  GTO  23 

653  X=V? 

703  RCL  32 

753*LBL  80 

654  GTO  61 

704  GTO  26 

754  38 

655  XEQ  60 

?85»LBL  63 

755  GTO  23 

656  teW 

786  RCL  33 

?56*LBl  81 

657  GTO  62 

78?  GTO  26 

757  35 

658  XEQ  80 

?08*LBL  64 

758  GTO  23  ■ 

£59  v= v? 

789  RCL  34 

759H6L  82 

668  GTO  63 

710  GTO  26 

768  98 

661  XEQ  00 

711*LBL  65 

761  GTO  23 

662  X=Y? 

712  RCL  35 

762*LBL  33 

663  GTO  64 

713  GTO  26 

763  95 

664  GTO  65 

714*lBL  66 

764  GTO  23 

665  GTO  26 

715  18 

765*LBL  84 

666*LBL  50 

716  GTO  23 

766  188 

66?  RCL  28 

71?*LBL  67 

76?  GTO  23 

668  GTO  26 

718  15 

76B»LBL  -pouhh* 

669»LBL  51 

719  GTO  23 

769  FIX  8 

678  RCL  21 

728HBL  63 

778  RhB 

671  GTO  26 

721  25 

771  FIX  4 

672»LBL  52 

722  GTO  2? 

772  RTN 

673  RCL  22 

?23*LBL  69 

773*LBL  * FIXE YEP 

674  GTO  26 

724  38 

774  STO  05 

675*LBL  53 

725  GTO  23 

775  2 

676  RCL  23 

726*LBL  70 

776  HOD 

677  GTO  26 

727  35 

777 

678HBL  54 

728  GTO  23 

7?3  RTN 

679  RCL  24 

729»LBL  71 

779  RCL  05 

638  GTO  26 

730  48 

788  1 

681HBL  55 

731  GTO  23 

781  * 

682  RCL  25 

732M.BL  72 

782  STO  05 

683  GTO  26 

733  45 

783  PTN 

684*L8L  56 

714  GTO  23 

734  .END, 

685  RCL  26 

735»L6L  73 

686  GTO  26 

736  58 

687*LBL  57 

73?  GTO  23 

688  RCL  2? 

738*LBL  74 

689  GTO  26 

739  55 

698+LBL  58 

740  GTO  23 

691  RCL  28 

74ULBL  75 

692  GTO  26 

742  68 

693*L8L  59 

743  GTO  23 

694  RCL  29 

744»LBL  76 

695  GTO  26 

745  65 

696»LBL  68 

746  GTO  23 

697  RCL  38 

747HBL  78 

698  GTO  26 

748  78 

699*LBL  61 

749  GTO  23 

780  RCL  31 

758»LBL  79 
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IV.  Contrail  Formation 

A.  Objective 

To  determine  whether  or  not  contrails  will  form  at  a  certain  level  in  the 
atmosphere  based  on  the  curves  from  the  Contrail  Forecasting  Manual  (Chief  of 
Naval  Operations,  1964). 

B.  Principles 

The  critical  temperature  determines  if  contrails  will  form  at  a  level. 
A  legist  squares  approximation  is  given  for  the  critical  temperature: 

Ter  it  =  ai  +  a2  In  p  +  a3  (In  p)2  +  a4  RH  +  as  (RH)2; 

f 

where: 

p  is  the  pressure  level; 

RH  is  the  relative  humidity; 

ai  =  -90.4994;  a2  =  3.4232;  aa  =  0.5587;  at  =  -0.0372;  as  =  0.0012. 

\ 

This  equation  was  determined  from  the  curves  in  Figure  9  of  the  Contrails 
Forecasting  Manual  (Chief  of  Naval  Operations,  1964). 

If  the  ambient  temperature  is  less  than  the  critical  temperature,  then 
contrails  will  form.  If  the  ambient  temperature  is  greater  than  the  critical 
temperature,  then  contrails  will  not  form. 

There  is  a  narrow  range  of  temperatures  where  the  relative  humidity  is 
important.  If  the  dew  point  is  known,  the  relative  humidity  can  be  calculated 
by  using  the  definitions  of  ..saturation  vapor  pressure  (e»)  and  relative 
humidity  (R.H.). 

e« (T)  =  6.11  exp  (L/Rw  (1/273  -  1/T)) 

R.H.  =  100  *  (e/e,) 

where  e*(T)  is  the  saturation  vapor  pressure 
e  is  the  actual  vapor  pressure 


T  is  the  temperature  (°K) 

L  is  the  latent  heat  of  evaporation  (J/kg) 

Rw  is  the  specific  gas  constant  for  pure  water  vapor  (J/kg  °K) 

If  the  dew-point  temperature  is  unknown,  the  relative  humidity  is 
estimated  and  a  +2°C  error  margin  is  assumed  when  computing  the  critical 
temperature.  If  the  ambient  temperature  is  within  +2°C  of  the  critical 
temperature,  a  "probably"  will  precede  the  "contrails"  or  "no  contrails" 
message.  The  following  criteria  are  used  to  estimate  the  relative  humidity: 

Taab  <  Tcrit  — >  (RH  =  0)  — >  contrails 
Taab  <  Tcrit  — >  (RH  =  100)  — >  no  contrails 
p  <  225  — >  (stratosphere)  — >  RH  -  0% 
p  >  300  — >  (not  in  upper  troposphere)  — >  RH  -  40% 

225  <  p  <  300  — >  (upper  troposphere)  — >  RH  —  40%  UNLESS 
cirrus  are  at  this  level  — >  RH  —  60% 
flow  is  known  and  from 

moist  region  — >  RH  —  60% 
dry  region  — >  RH  -  0% 


C.  User  s  Guide 

1)  Make  sure  both  the  calculator  and  printer  are  on  and  that  the  printer  is 
set  to  the  MAN  mode. 

2)  |  XBQ 1  1  alpha!  SIZE  1  alpha 

3)  Prompts  SIZE _ ;  Type  100. 

4)  Load  all  three  cards,  sides  1-6. 

5)  After  "WORKING”  signal,  program  is  ready  for  execution. 

6)  f  xEQ  1  |  alpha  1  TRAILS  I  alpha 

7)  Program  asks  about  CHANGING  FLAGS.  For  the  first  run  through,  leave  the 
flags  as  they  are  and  continue  by  hitting  the  1  R/S  |  key. 

8)  Enter  a  pressure  level  in  mb. 

9)  Press  |  R/S|  key.  (Hit  | R/S |  after  each  entry. ) 

10)  Enter  a  temperature  in  °C. 

11)  A  prompt  for  a  dew-point  temperature  will  be  displayed.  If  Td  is 

known,  type  a  1;  if  Td  is  unknown,  type  a  0. 

12)  The  program  will  then  decide  whether  or  not  contrails  can  be  formed  at 

that  particular  level  depending  upon  the  criteria  discussed  on  the 

preceding  pages.  If  the  dew-point  temperature  is  unknown,  there  might  be 
more  "yes"  or  "no"  questions  depending  upon  the  pressure  level  and 
ambient  temperature.  Answer  "1"  for  "yes"  and  "0"  for  "no." 

13)  The  program  will  determine  whether  or  not  contrails  .* will  form  at  a 

certain  pre-established  pressure  level. 

14)  After  the  answer  is  displayed,  execution  can  continue  for  another  level 
by  pressing  |  R/S] 


By  changing  the  flags,  two  options  are  possible.  Setting  Flag  1  will 
make  the  prompts  shorter  and  less  specific.  Setting  Flag  2  means  that 
the  program  assumes  that  the  dew-point  temperature  is  known  and  skips  the 
questions  concerning  relative  humidity  estimation. 

To  change  a  flag:  a)  Press  gold  button  and  [ SF  |  key. 

b)  Program  prompts  SF  . 

c)  Enter  either  (or  both):  01  -  short  prompt 

02  -  known  Ta 

d)  To  clear  flags,  press  gold  button  and  [  CF  |  key. 

e)  Program  prompts  CF _ . 

f)  Enter  01  or  02  to  clear  any  flag  that  was  set. 

The  program  assumes  all  flags  are  clear  each  time  it  is  first  executed. 

D.  Examples 

Examples  A-E  in  Table  IV  illustrate  the  use  of  the  program  "TRAILS"  in 
deciding  whether  contrails  will  form  from  upper-level  sounding  data.  In 
example  A,  where  no  flags  have  been  set,  a  pressure  of  450  mb  and 
temperature  of  -35° C  are  entered  in  the  "TRAILS"  program.  The  dew-point  is 
unknown  in  this  case  so  a  "0"  is  entered.  The  program  decides  that  contrails 
will  not  be  possible  with  these  conditions  which  agrees  with  the  contrail 
curve  diagram.  (See  Figure  9.)  Example  B  is  similar  to  A  (no  flags  set  and 
dew  point  unknown),  but  contrails  will  probably  form  with  the  particular 
upper-level  sounding  data.  In  example  C,  no  flags  are  set,  but  the  dew  point 
is  known.  With  a  pressure  of  250  mb,  a  temperature  of  -60.0°C,  and  dew  point 
of  -62.0°C,  contrails  will  form.  Examples  0  and  E  have  flags  set,  which 
result  in  shortened  prompts  for  both  cases,  and  for  example  E,  the  dew  point 
is  assumed  to  be  known.  Both  examples  0  and  E  illustrate  that  contrails 
will  not  form  in  these  environmental  conditions. 


Table  IV 

Examples  of  the  Contrail  Formation  Program. 


A  B 

p=’  m 


P=’  H8 

P=?  mb 

250.0000 

*** 

450.0000 

»** 

T=’  C 

T=?  C 

-54.0000 

*** 

-35.0000 

**» 

DO  YQU  KNOW 

BO  YQU  KNOW 

DEWPOINT,  TD’ 

DEWPOINT-  TB? 

YES=0NE 

YES=0NE 

N0=ZEP0 

H0=ZER0 

KNOW  TD’ 

KNOW  TD? 

0.0000 

*** 

0.0000 

*** 

ARE  THERE 

NO  CONTRAILS 

CIRRUS  AT 

THIS  LEVEL’ 

No  flags  set 

Td  unknown 

YES=0NE 

N0=ZER0 

CIRRUS’ 

1.0000 

*** 

PROBABLY 

CONTRAILS 

No  flags  set,  Td  unknown 

C  D 


P=’  HB 

P=?  H8 

300.0008 

*** 

225.0000 

*** 

T=’  C 

T=’  C 

■■ 

-48.0000 

*** 

-60.0000 

*** 

KNOW  TD? 

DO  YOU  KNOW 

0.0008 

*** 

DEWPOINT,  TD’ 
YES=ONE 

CIRRUS’ 

8.0000 

*** 

NO=ZERO 

KNOW  TD’ 

KNOW  FLOW’ 

1.0000 

*** 

1.0000 

*** 

HST  OR  DRY’ 

TD=’  C 

1.0000 

*** 

-02.0000 

*** 

NO  CONTRAILS 

CONTRAILS 

No  flags  set 
T  d  known 


Flag  01  set 
Unknown  T  d 


P=?  H8 

850.0000 

*** 

T=’  C 

-10.0000 

Mt 

TD=?  C 

-20.0000 

*** 

HO  CONTRAILS 

Flag  02  set 
T  d known 


1000 


rj  -  .  sr_  «rj  vucj 


»i 


I  i 

»J 


Program 


0!*L8l  -TRAILS’ 

82  CF  hi 

93  CF  02 

94  ‘CHANGE  FLAGS’1 
05  PROMPT 
06«L8L  10 

07  TONE  l 
08  -P=>  MB’ 

0?  FS?  21 

10  PR  A 

11  PROMPT 

12  FS''  21 

13  PRX 

14  STO  80 

15  IN 

16  STO  01 
1?  TONE  2 

18  *T-’  C* 

19  FS?  21 
28  PRA 

21  PROMPT 

22  FS?  21 

23  PRX 

24  STO  04 

25  FC?  62 

26  GTO  03 
27*L8L  01 

28  XEQ  -RH‘ 

29  XEQ  *TC‘ 

30*LBL  02 

31  RCL  04 

32  RCL  83 

33  X’Y’ 

34  GTO  C 
35*LBL  B 

36  "NO  CONTRAILS* 

37  TONE  5 

38  AVI  EH 

39  ADV 

40  STOP 

41  GTO  10 
42*LBL  C 

43  ‘CONTRAILS’ 

44  TONE  8 

45  AVI EH 

46  ADV 

47  STOP 

48  GTO  10 
49»L8L  03 
50  FS?  81 


Table  V  Contrail  Formation  Program. 


51  GTO  11 

101  STO  82 

151  XROH  28,33 

52  -DO  YQU  KNOW 

182  FS’  01 

152  X=0? 

53  AVI EH 

103  GTO  12 

153  GTO  A 

54  PSE 

184  -ARE  THERE- 

154  FS’  01 

55  -DEHPOINT,  TP?- 

185  AVIEH 

155  GTO  14 

56  AVIEH 

1@6  PSE 

156  -FROH  HOIST- 

57  PSE 

107  -CIRRUS  AT- 

157  AVIEH 

58  ‘YE;  E‘ 

108  AVIEH 

158  PSE 

59  AVIE* 

109  PSE 

159  -REGION- 

60  PSE 

110  -THIS  LEVEL’* 

160  AVIEH 

61  -N0=2ER0- 

111  AVIEH 

161  PSE 

62  AVIEH 

112  PSE 

162  •=  ONE- 

63  PSE 

113  -YES=OME- 

163  AVIEH 

64»LBL  11 

114  AVIEH 

164  PSE 

65  -KNOH  TD?‘ 

115  PSE 

165  'FROH  DRV 

66  FS?  21 

116  *NO=ZEPO" 

166  AVIEH 

67  PRA 

117  AVIEH 

167  PSE 

68  XROH  28,21 

118  PSE 

168  *=  ZERO- 

69  XROH  28,33 

119*L8L  12 

169  AVIEH 

78  CHS 

120  -CIRRUS’* 

170  PSE 

71  1 

121  FS?  21 

171*LBL  14 

72  X=Y? 

122  PRA 

172  -HST  OR  DRY’ 

73  GTO  81 

123  FS?  21 

173  FS?  21 

74  8 

124  XROH  29,33 

174  PRA 

75  STO  82 

125  PRX 

175  XROH  28,33 

76  XEQ  *TC* 

126  X*Y? 

176  5 

77  RCL  04 

127  GTO  A 

177  PRX 

78  RCL  03 

128  XROH  28,33 

178  * 

79  X>Y’ 

129  48 

179  GTO  39 

38  GTO  C 

130  STO  02 

180  RCL  01 

81  100 

131  FS?  01 

181  X<9? 

82  STO  82 

132  GTO  13 

182*LBL  P 

83  XEQ  -TC- 

133  -DO  YOU- 

183  XEQ  -TC* 

84  RCL  03 

134  AVIEH 

184  PCL  03 

85  RCL  04 

135  PSE 

185  RCL  04 

86  X>Y’ 

136  -KNOH  FLOW?’ 

186  - 

87  GTO  B 

137  AVIEH 

187  ABS 

88  8 

138  PSE 

188  2 

89  STO  82 

139  -YES=ONE- 

189  X<=Y? 

98  RCL  00 

140  AVIEH 

198  GTO  82 

91  225 

141  PSE 

191  -PROBABLY- 

92  X>Y? 

142  -NO=ZERO* 

192  AVIEH 

93  GTO  A 

143  AVIEH 

193  PSE 

94  48 

144  P$E 

194  GTO  02 

95  STO  82 

145*LBL  13 

I95*LBL  -RH- 

96  380 

146  ‘KNOH  PLOW?- 

196  RCL  04 

97  RCL  08 

147  FS?  21 

197  273 

98  X>Y? 

148  PRA 

198  ♦ 

99  GTO  A 

149  FS’  21 

199  XEQ  -ES- 

108  68 

150  PRX 

288  $T0  02 

201  tQnE  T 
282  ’TD=?  C- 

203  FS’  21 

204  PRA 

205  FS?  2! 

206  PRX 

287  -XR0K  28,33 

288  273 

209  * 

210  XEQ  -ES- 

211  RCL  82 

212  •' 

213  180 

214  * 

215  $T0  02 

216  RTN 
217»LBL  -ES* 

218  1/X 

219  CHS 

228  273 

221  1/X 

222  * 

223  2580 

224  * 

225  8.461 

226  / 

227  EtX 

228  6 

229  * 

238  RTN 

231 ♦LSI  -TC- 

232  -98.4994 

233  STO  83 

234  .558? 

235  RCL  01 

236  * 

237  3.4232 

238  ♦ 

239  RCL  01 

248  * 

241  ST+  03 

242  .0012 

243  RCL  02 

244  * 

245  -.8372 

246  ♦ 

247  RCL  82 

248  * 

249  ST*  03 

250  RTN 

251  .END. 


*4l 


i 


WA  V»  , 


V. 


Determination  of  Lifting  Condensation  Level 
Objective 

Given  a  surface  temperature,  surface  pressure  and  specific  humidity,  this 


A. 


program  calculates  the  lifting  condensation  level. 

B.  Principles 

The  lifting  condensation  level  (LCL)  is  the  level  to  which  a  parcel  of 
air  can  be  lifted  dry  adiabatically  before  it  becomes  saturated.  During  the 
lifting  process,  the  potential  temperature  of  the  air  parcel  and  the 
saturation  mixing  ratio  remain  constant.  The  actual  mixing  ratio  decreases 
and  eventually  equals  the  saturation  mixing  ratio.  The  lifting  condensation 
level  is  obtained  when  the  actual  vapor  pressure  equals  the  saturation  vapor 
pressure  (e  =  e» ) .  Integrating  the  Clausius-Clapeyron  equation  from  273  K  to 
air  temperature  (T) 


y 

273 


den 

es 


T 

= ; 

273 


L 

Rw 


dT 

T2 


where  e«  is  the  saturation  vapor  pressure;  L  is  the  latent  heat  of 
condensation  (2.5  x  10®  J/kg) ;  Rw  is  the  gas  constant  for  water  vapor  (462 
J/kg  °K);  and  T  is  the  temperature  (°K),  results  in  the  following  expression 
for  the  saturation  vapor  pressure: 


-  =  6-u  exp  (sr  '^3  -  i»- 

Combining  the  hydrostatic  equation  and  the  ideal  gas  law  and  expressing  the 
temperature  as  To  -  Y  z  results  in  the  following  expression  relating  pressure 
(P)  to  height  (z): 


where  y  is  the  dry  adiabatic  lapse  rate,  To  is  the  standard  atmospheric 
temperature,  w  is  the  specific  humidity  and  Po  is  the  standard  atmospheric 
pressure,  (w  -  .622  (e/P)  ==>  P  =  .622  e/w) 


s 


An  expression  for  the  actual  vapor  pressure  (e)  is  obtained  by  using  the 
definition  of  specific  humidity  and  substituting  the  above  expression  in  for 
P.  The  lifting  condensation  level  is  obtained  when  the  following  two 
equations  are  equal  to  each  other: 


-Y 

wPo  (To  -yz  )R 
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After  a  series  of  substitutions,  using  Taylor  series  expansion  of 
In  x  ~  x— 1  and  several  algebraic  manipulations,  the  following  equation  is 
derived  giving  the  lifting  condensation  level: 

rjpQ  wPq 

=  T.  (102.041  -  14.6429  (23.3058  -  -  [l^SoS5 


-  23.30582  -  V  )) 

To 


Only  surface  data  values  of  temperature  (To),  pressure  (Po)  and  specific 
humidity  (w)  are  required  to  find  the  LCL. 


C.  User’s  Guide 

1)  Turn  on  calculator  and  set  printer  to  the  NOR  mode. 

2)  [  XEQ|  LCL 

3)  Prompt  reads,  SFC  TEMP  °C? 

4)  Type  in  a  surface  temperature. 

5)  Press  R/S  to  enter  the  temperature. 

6)  Prompt  is  SFC  PRES  MB? 

7)  Enter  a  surface  pressure. 


9)  Prompt  reads,  SP  HUMID  G/G? 

10)  Enter  a  specific  humidity  value. 

11)  1  R/S 

The  lifting  condensation  level  will  appear  with  its  value  in  meters.  If 
the  LCL  occurs  at  the  surface,  the  word  SURFACE  will  be  displayed. 

D.  Examples 

Examples  A,  B  and  C  in  Table  VI  illustrate  the  use  of  the  program  "LCL" 
in  estimating  the  lifting  condensation  level  from  a  surface  temperature, 
pressure  and  specific  humidity.  In  example  A,  a  surface  temperature  of  9°C, 
a  surface  pressure  of  1010  mb,  and  a  specific  humidity  of  .006  g/g  are  entered 
in  the  "LCL”  program.  The  data  result  in  a  lifting  condensation  level  of 
318.9m  which  is  verified  by  a  thermodynamic  diagram.  (See  Figure  10.)  In 
example  B,  the  lifting  condensation  level  occurs  at  the  surface  and  in  example 
C,  the  LCL  is  at  637.4  m. 
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Table  VI 

Examples  of  the  Lifting  Condensation  Level  Program 


B 


$*C  TEH?  C? 

9.0000 

SrC  PRES  MB'7 


SFC  TEMP  iy 

.5000 

SFC  PRES  HB7 


1/010.0000 
SP  HUMID  S/G? 

.  0060 

LCL  H. =3 IS, 895? 


1  ■008.8808 
SP  HUH  ID  Q/ip 

.0040 

SURFACE 


C 


SFC  TEMP  rj 

6.0888 
SFC  PRES  HE? 

1/020.0000 
SP  HUMID  G/G? 

.0040 

LCL,  H. =637.3732 
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Figure  10 

Lifting  condensation  levels  determined  by  surface  temperature,  pressure 
and  specific  humidity. 
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E.  Program 
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Table  VII 

Lifting  Condensation  Level  Program. 


91 *181  "LCL* 

82  23.3058 
03  '-TO  25 
84  -SPC  TEMP  V- 
05  PROMPT 


06  273.15 
87  ♦ 

8?  STO  2! 

80  -SFC  PRES  MB'7* 
18  PROMPT 


12  *SP  MIJMI5  r, 
1?  PROMPT 

14  STO  23 

15  RCL  22 

16  * 

17  3.38842 
IS  ' 

19  LM 
23  STO  24 
21  RCL  25 


23  ?;+2 

24  75418.2 

25  RCL  2’. 

26  / 


28  SORT 


28  STO  26 

30  CMS 

31  RCL  24 


34  RCL  25 


36  14.6429 


182.041 


4  j  PCL 


42  * 

43  STO  38 

44  0 

45  XOV 

46  X<=V> 

47  QTO  -SRC* 

48  "LCL..  M,.=* 

49  fiRCL  30 

50  MV I £M 

51  GTO  *EH8* 
52*LBL  “SFC* 

53  "SUP^QCt* 

54  fiVJEH 
55+LBL  "EMD* 
56  .END. 
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